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The formal link between the linear combination of atomic orbitals approach to 
density functional theory and two-center Slater-Koster tight-binding models is used 
to derive an orthogonal d-band tight-binding model for iron with only two fitting pa- 
rameters. The resulting tight-binding model correctly predicts the energetic ordering 
of the low energy iron-phases, including the ferromagnetic BCC, antiferromagnetic 
FCC, HCP and topologically close-packed structures. The energetics of test struc- 
tures that were not included in the fit are equally well reproduced as those included, 
thus demonstrating the transferability of the model. The simple model also gives 
a good description of the vacancy formation energy in the nonmagnetic FCC and 
ferromagnetic BCC iron lattices. 

PACS numbers: 71.20.Be,75.50.Bb,71.15.Ap 

I. INTRODUCTION 

While Kohn-Sham (KS) density functional theory (DFT)^ has found very broad appli- 
cation for the simulation of interatomic bonding, its computational cost still places limita- 
tions in its application when treating the length scales necessary for the strain fields from 
dislocations^ or light elements in metals.- Furthermore as the system size grows the number 
of configurations needed for thermodynamic integration becomes intractable. This makes 
the use of computationally efficient parameterized methods attractive. The continued inter- 
est in parameterized methods also comes from the obvious wish to gain physical insight. In 
this respect one of the most successful methods is the tight-binding (TB) method. 

In its conventional form the TB method models the total energy as a repulsive pair 
potential and a bonding many-body term. The bonding energy is obtained by solving a 
two-center Slater-Koster (SK) Hamiltonian.- Following TBs empirical introduction several 
conceptual advances, mainly the TB bond model^^S^ the Harris-Foulkes functionaF'^ and 
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the related second-order expansion of the KS energy^"— have been made. Together these 
provide an appeahng conceptual framework, but in practice there are several "philosophies" 
on how the parameterization should be performed and the success of TB depends on this 
parameterization.—"— 

There is thus a demand for TB parameterizations based as closely as possible on the 
DFT energy functional. In the present paper we construct an orthogonal TB model for iron. 
Special focus is put on using a limited number of fitting parameters without compromising 
the predictive quality of the model. We demonstrate how the formal link between DFT 
linear combination of atomic orbitals (LCAO) methods and two-center TB may be used to 
obtain the TB bonding energy. This is achieved by down-folding a pseudo-atomic orbital 
(PAO) basis onto a minimal basis set. We demonstrate the transferability of both basis 
functions and bond-integrals, thereby validating the two-center approximation. We show 
how the resulting TB model for iron correctly predicts the energetic ordering of the low 
energy iron-phases, including the ferromagnetic (FM) BCC, antiferromagnetic (AFM) FCC 
and topologically close-packed structures. Finally, we test the transferability of the model 
on the vacancy formation energy in the NM-FCC and FM-BCC iron lattices. 



II. METHOD 
A. Background 

In LCAO the basis functions are written as product of a radial part with an angular 
function 

= 0/iim(r) = Uiji{r)Yim{r) (1) 

We use the capital indexes / and J to label atoms and the index as a condensed index for 
the angular character Im. We leave out the principal quantum number as we only treat the 
valence states. While minimal basis sets use just one basis function for each valence atomic 
orbital, the variational flexibility of LCAO basis sets can be improved by adding several 
radial functions for a given angular momentum, so-called multiple-^ basis functions. The 
index j in Eq. ([T]) counts the number of radial functions for a given angular character /i. 
Furthermore higher spherical harmonics, so-called polarization functions, are often added 
to further improve the basis. By expanding the Kohn-Sham (KS) orbital wave functions in 
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FIG. 1: Illustration of the down- folding of a triple-C basis to an optimal single-^ basis. The 
original 3 — C GPAW pseudo-atomic orbitals (PAO) basis is shown to the left. The plot to the right 
shows the optimal basis function for Fe in the simple cubic structure (with a lattice constant of 
a = 2.50 A), the FCC (a = 3.46 A) and the BCC (a = 2.87 A) structures. The structures all have 
a nearest neighbour distance of 2.5 A and the basis functions are virtually indistinguishable. The 
confinement potentials corresponding to AEp^o = 0.1 eV are shown in black. Also shown with a 
dashed line is the optimal basis function for the Fe dimer at an interatomic distance of 2.5 A. 

terms of a basis set 

l^n) = E4"ilfe) (2) 

the KS equations can be written in matrix form, which introduces the Hamilton and overlap 
matrices 

Jju Jju 

In the present paper we will use the radially confined PAOs^ implemented in the GPAW 
code for the radial functions in Eq. ([1]).— The PAO basis functions have a well defined 
radial extent due to the confinement potential used, see Fig. [H— "i^. Confining the radial 
extent of the atomic orbitals increases their energy. Following the original work— this 
energy shift, /S.EpAOi is used to define the radial cut-off. For most part of the paper we use 
the standard setup of GPAW, AEpAo = 0.1 eV, which leads to confinement radii of 4.7 A 
for the s-PAQ and 2.7 A for the rf-PAO of iron, and an onset of the confining potential at 
60 % of the confinement radius. 
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FIG. 2: Comparison of the density of states of non-magnetic iron calculated using three different 
basis sets. The lattice constants for the calculations were FCC a = 3.46 A and BCC a = 2.87 A. 
The structures have a n.n. distance of 2.5 A. 

In order to achieve the precision of a systematic grid or plane wave basis, an atomic 
basis must include both multiple-(^ and polarization basis functions, thus far removed from 
the simple TB models that we wish to construct. We therefore use the dual basis sets of 
grid points^^ and atomic orbitala^^ implemented in the GPAW code. We first calculate self- 
consistent total energies and potentials using the systematic grid basis. We then obtain the 
eigenstates {ipn) expanded in a 3-^ basis, Eq. (|2]), by performing a single diagonalization in 
the potential obtained by the grid calculation. Fig. [2] illustrates the very good agreement 
between the DOS calculated with the grid basis and with a 3-^ basis. 



B. Optimized Atomic Orbitals 

The optimized minimal (1-C) basis is obtained from the multiple-^ basis by a down-folding 
of the LCAO eigenstates for a given atomic configuration. In a non-orthogonal minimal basis 
{\(pi^)}, the contravariant basis {(v^'*!} provides a simple expression for the closure relation 

{^"'\ = J2s7;jM , $^i¥./.)(/i = i (4) 

Ju IfJ, 

with the overlap matrix S = {(fjf^\(fj^). The closure relation may be seen as a projection 
operator, which if applied on | ■?/'„), measures to which extent | can be represented in 
the basis. We thus write the projection of \ipn) expanded in the multiple-^ basis 
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Eq. ([2]), on the minimal basis {{fi^)} as 

Pn = J2{^n\^I,){^"'\^n) , P = iV"! /„P„ (5) 

Ifi n 

where fn is the occupation of the eigenstate n and Ne the number of valence electrons. The 
basis function is written as a linear combination of the 3-^ basis-functions for the same 
angular character 

= ^Oiiji(f)jjf,{r) (6) 
j 

The coefficients ajji, Eq. ([6]), are found by maximizing the projection P, Eq. ([5]). Eq. ([5]) 
was introduced earlier for reducing multiple-^— and plane wave basis sets^^ to minimal basis 
sets. It has however not been broadly applied for this purpose because the optimal basis 
for a given structure is not transferable. This is less of a problem for TB where we wish 
to parameterize the bond integrals as a function of interatomic distance. Fig. [1] shows that 
for a given interatomic distance there is a very good agreement for the 3(i-PA0 between the 
two extreme cases of a close packed solid Fe and the Fe dimer. For the 4s-PA0 there is also 
a very good agreement between the solids, whereas the 4s orbital for the dimer contracts 
somewhat. 

Eq. ([5]) was first used for defining optimal AOs for TB from a plane wave basis by Meyer 
and coworkers.—"— Our method differs through the choice of an LCAO basis for which 
makes the down-folding a numerical simpler procedure. Eq. can be calculated using only 
the variational coefficients c''"\ the overlap matrix and the sparse matrices containing the 
coefficients a, Eq. ([6]). We maximize P with respect to a using a standard conjugate gradient 
method and have found the same minimum for all test cases irrespective of starting values. 
A further feature of the present method is that the basis underlying the TB parameters has 
a well defined radial extent meaning that its infiuence on the bond integrals may be studied 
systematically. 

Constructing a minimal sd-basis for the FCC and BCC-iron structures used for Fig. [2] 
gave P = 0.995 for both. Not surprisingly P ^ 1 also means that the DOS calculated with 
an optimized basis is very similar to the 3-( DOS. We have also compared to the DOS found 
by optimizing the band energy directly and found it virtually indistinguishable from that 
obtained through projection. 
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C. TB Energy Functional 

To a good approximation the structural energy of the transition metals is determined by 
the d-valence^ while the contribution of the s-electrons may be approximated by a volume 
dependent embedding contribution. For the evaluation of the TB energy we further assume 
that the charge transfer in Fe is small and may be neglected. We therefore assume that the 
atoms remain charge neutral and only allow for magnetic fluctuations, such that our TB 
energy functional is given as 

ExB ^bond ~l~ ^mag ~l~ E^ep ~\~ E^^nb E j^ee— atoms (7) 

The first term is the bond energy of the d-electrons within the TB bond model^i^ which for 
coUinear spins may be written a^ 

Ebond = X] X] P'lfiJi'Hlf.Ju (8) 

where a labels the spin. As we assume local charge neutrality the second-order term of 
the expansion of the DFT energy only contains a magnetic contribution depending on the 
Stoner exchange integral.— The second term in Eq. ([7]) is the Stoner exchange energy^Si"— 

Emag = -^^Ijm^J (9) 

J 

where mj is the magnetic moment on atom J . We further approximate the Stoner parameter 
I J as an atomic quantity. The third term in Eq. ([7j) is a pair-wise repulsive contribution 
modelling the double counting term of the TB bond energy.- We write the repulsive potential 
simple exponential 

Erep = <ip^M-KipRij) (10) 

I,J^I 

Finally, Eq. ([7]) approximates the contribution of the s-electrons to the cohesive energy with 
a simple embedding term. Based on the second-moment approximation to the DOS, we 
model this as having a square-root dependence on the coordination number, n = 1/2.—"— 



E, 



= - E f ^M-KLRh)) (11) 

n = 1 would correspond to a pair potential. For the embedding function we use a Gaussian 
like radial dependence. This has been proposed earlier— and will be justified later in this 
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paper. Finally, the term Efree-atoms corresponds to the energy of the atoms at infinite 
separation. 

D. Bond Integrals 

We have calculated the band structure for a series of interatomic distances for the iron 
dimer and for iron in the FCC and BCC structures. The calculations were performed by first 
calculating a self-consistent potential using the grid basis of GPAW.^° Then a diagonalization 
was performed using a standard 3-^ PAO basis of the GPAW— which was then down-folded 
in a minimal basis by maximizing the projection, Eq. (jS]). 

For a S(i- minimal basis 6x6 sub-matrices of the LCAO Hamilton Hi^jy or overlap Sj^^jy 
matrices are associated with each pair of atoms. Each of these matrices can be rotated into 
a bond-oriented coordinate system, resulting in the bond-integrals 

Pif^ju = '^U^f.jy'Hiyij^iUi^'ju (12) 

where Ui^j^ is the matrix that rotates the global coordinate system into a bond-oriented. 
In the two-center approximation,^ by symmetry only the sscr, sda, ddcr, ddir and dd6 matrix 
elements are non-zero. In our orthogonal rf-valent TB model we will retain only the dda, 
ddiT and dd6 integrals. 

In Fig. |3] we show the bond-integrals f3 that were calculated from the optimal minimal 
basis using Eq. f lT2|) . The bond integrals are discontinuous and poorly transferable. It has 
earlier been shown that including screening makes the bond-integrals /3 continuous at the 
n.n. and n.n.n. distances. -^i^i^i^^ This prompted us to define the bond-integrals based on a 
Hamiltonian orthogonalized by a symmetric Lowdin procedure,— 

H = S-^/^HS-^'^ (13) 

where H corresponds to the full Hamiltonian in the sd minimal basis. Compared to other 
orthogonalization schemes the Lowdin orthogonalization has two important advantages: the 
orthogonal orbitals bear the same symmetry as the non-orthogonal original vectors,- and 
are the closest in a least squares sense.— Fig. |3b shows that the bond-integrals obtained by 
using H in Eq. (fT2|) are both transferable and continuous. The very good agreement shown 
in Fig. |3]d even with the Fe-dimer is somewhat surprising. It has already been shown in 
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Fig. [T] that the optimal d-basis is transferable for a given interatomic distance. Therefore 
the poor transferability observed in Fig. [3^ can only be due to three-center, l^ly^j), 
contributions to the Hamilton matrix elements leading to an environmental dependence of 
the two-center integrals. The effect of the Lowdin orthogonalization must be a screening of 
the three-center integrals. 

A qualitative rationalization of the transferability can be found by comparing H to the 
D matrix used in an analysis of chemical pseudopotential theory. Large three-center con- 
tributions will be associated with large two-center overlap integrals thereby screening the 
large three-center integrals. This interpretation is confirmed in Fig. |3fc where radial extents 
of the basis functions, and thereby the three-center contributions, are reduced. Using a 
AEpAo = 0.4 eV instead of AEpAo = 0.1 eV reduces the radial extent of the (i-orbitals 
from 5.1 A to 3.9 A. Consequently the unscreened bond-integrals show transferability and 
are continuous. 

The bond integrals are fitted to simple exponentials as 



Due to the transferability of the bond-integrals. Fig. [3], we simply use the bond-integrals 
obtained for the dimer, the parameters are given in Table [H At the nearest-neighbour 



integrals dda : ddn : dd6 = — 0.60eV : 0.41eV : — O.OSeV shows a surprisingly good agreement 
with the canonical d-hand ratio of —6 : 4 : —1.— The transferability to the dimer also forms 
a link to the widely used DFTB approach^, where the bond-integrals are evaluated from a 
dimer calculation using a single-^ basis in a potential from overlapping atomic densities.— 
To a certain degree Fig. [3] may be seen as a validation of this approach. However, it should 
be pointed out that the transferability obtained in Fig. [3b holds only for the short-ranged 
(i-orbitals. The longer-ranged s-orbitals will be the subject of a future study. To this end 
the fact that our matrix elements are evaluated in the actual crystal potential is a clear 
advantage when studying the influence of three-center integrals. 
A cut-off function given as 



Pdd\{.R) = addxexp{-bMxR) 



A = 0", TT, (5 



(14) 



distance of the BCC and FCC structure of around 2.5 A the relative strength of the bond 



f{R) = I \ (cos (7r( 



R—[Rcut—dcut) 
dcut 



1 



)) +i) .R, 




'Cut 



(15) 







, R ^ R, 



■cut 



9 




R(A) 



FIG. 3: Bond integrals: a) non-orthogonal AEpAo = 0.1 eV. b) orthogonal AEpAo = 0.1 eV. 
c) non-orthogonal AEp^o = 0.4 eV. The full lines in b) show a fit to simple exponentials to the 
orthogonal AEpAo = 0.1 eV dimer curves. 
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dda 




-34.811 


1.625 




ddiT 




63.512 


2.014 




ddS 




-50.625 


2.597 


dcut 
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(A) 


0.5 


3.5 




E 




1031 


3.25 




Eemb 




3.70 


0.23 


dcut 


Rent 


(A) 


0.5 


5.5 



TABLE I: Parameters of the tight binding model, Eqs. ([TO]) . (fTT]l and (fHll . The units of hemh are 
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FIG. 4: Density of states obtained with the orthogonal d-band model. The structures are as in 
Fig. El 



was applied to the distance-dependent pair-interactions. The cut-off parameters are given 
in Table |Tl] and were chosen so that the bond-integrals and pair and embedding potentials 
are cut-off around the onset of the d and s confining potentials respectively. The resulting 
DOS of the TB model are shown in Fig. HI Apart from the obviously lacking peaks due to 
sd-hybridization there is some disagreement with respect to the magnitude of the DOS at 
the Fermi-level. A good agreement is found between the location of the peaks. 

Omitting the s-electrons in the bond energy means that the number of d-electrons must 
be introduced as a parameter. As the FCC and HCP structures have the same first and 
second nearest neighbor shells, we assume that the embedding and repulsive energies for the 
two structures at equal volume is the same and the energy difference is purely due to the 
difference in Ebond- We thus use the energy difference of the FCC and HCP structures at 
equilibrium volume to fix A'^^ = 6.8 e/atom. Thereby a bond energy difference between the 
FCC and HCP structure of -53 meV in good agreement with the DFT value of -60 meV is 
obtained. 

Compared to earlier TB models of irort^"— i^"— our treatment of magnetism is similar to 
that of refs. y andE 



Instead of obtaining the Stoner exchange integral directly from DFT, 
we set it to / = 0.76 eV to get a good energy difference between the magnetic and non- 
magnetic structures. This choice leads magnetic moment of 2.65 /i^/atom and 1.34 /i^/atom 
at the equilibrium volumes of BCC Iron and FCC Iron respectively. Compared to DFT, 
2.21 /i^/atom and 1.05 /i^/atom, the magnetic moments found with our TB model are 
to large. We attribute this to the lack of S(i-hybridization in the model and see this as a 
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fundamental limitation of the present approach. Finally, we have tested the stability of the 
FM-BCC structure in our TB model by doing 500 MD steps at 300 K using a Andersen 
thermostat and a Velocity Verlet integrator. We find the FM-BCC structure to be stable. 

E. Repulsive and Embedding Energies. 

For the repulsive and embedding terms, Eqs. (TTOl)-( |TT|) . the exponents are fixed by the 
extracted bond and overlap integrals. The repulsive part we see as an overlap repulsion which 
should thus be proportional to the square of the most long-ranged dd-overlap integral. Using 
Pddcr = 1-625 A^-*^, suggest that we set brep = 3.25 A^^ The embedding part we see as arising 
from not including the s-states in the bonding term, it is thus written in terms of the square 
of the Pssa matrix element for the Fe2 dimer. We find this to be well represented by a 
Gaussian with an exponent of 0.115 which suggests 6emfe = 0.23 A"^. We thus end up 
with a TB-model where only two parameters must be found by fitting total energies. We 
fit the parameters a^ep and aemb, Eqs. (niJ]) -( lTT]) . to the DFT energy- volume curves for non- 
magnetic BCC, FCC and HCP structures. The resulting parameters are given in Table [H 
The resulting bulk moduli and phase stabilities are given in Table Ull Table Ull also shows 
the results of applying the TB-model to a number of topologically closed packed phases^ 
and the AFM-FCC and FM-BCC structures. It is seen that the agreement is similar to 
the structures included in the fit which demonstrates the transferability of the model. The 
main disagreement is the bulk modulus of the FM-BCC iron phase which is underestimated. 
We attribute this to the too large magnetic moment found with J = 0.76 eV leading to a 
high-spin state at extended volumes. 

F. Transferability 

We further test the transferability of the model by evaluating the vacancy formation 
energy (VFE) in FM-BCC and NM-FCC iron and the formation energy with respect to 
the solid of an NM-FCC-(lll) unsupported monolayer of Fe. The VFE are calculated in a 
2x2x2 cubic supercell, which thus holds 15 atoms for BCC and 31 for FCC. As shown 
in Table IIIII we find a reasonable agreement with DFT. In all three cases we find that the 
open structure is to low in energy, compared to the close packed. One would expect that an 
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increase in n in the embedding function, Eq. (fTT]) . would stabilize the close packed structure 
compared to the open. Consequently, we find that using an exponent of n = 0.55 instead of 
a square-root potential gives a better agreement with DFT for the formation energies of the 
open structure. Setting n = 0.55 and reoptimizing aemb and a^gp, again only fitting to the 
NM-BCC, NM-FCC and NM-HCP structures, we find a^ep = 1088 eV and aemp = 3.18 eV. 
The reoptimization can be done without changing the agreement found in Table HIl which 
shows that by introducing a more fiexible potential better agreement can be achieved at the 
expense of the simplicity of the model. 

III. CONCLUSION 

We have shown how to derive an orthogonal d-band TB model for iron with only two 
fitting parameters. The resulting TB model correctly predicts the energetic ordering of the 
low energy iron-phases, including the ferro-magnetic BCC, anti-ferromagnetic FCC and the 
topologically closed packed structures. We have found that test structures that were not 
included in the fit are equally well reproduced as those included, thus demonstrating the 
transferability of the model. The simple model gives a good description of the formation 
energy of a vacancy in the NM-FCC and FM-BCC iron lattices. 

Simple orthogonal TB models form the basis of the bond-order potentials (BOPs),— 
which in their simplest second-moment approximation are described by many-body energy 
terms that correspond to a square-root embedding function.— At the same time the BOPs 
constitute a systematic approximation of the TB model by including higher moment con- 
tributions to the binding energy. The present work could form a crucial link between DFT 
and interatomic potentials in a hierarchy of controllable accuracy. 

IV. ACKNOWLEDGMENTS 

We acknowledge financial support through ThyssenKrupp AG, Bayer MaterialScience 
AG, Salzgitter Mannesmann Forschung GmbH, Robert Bosch GmbH, Benteler Stahl/Rohr 
GmbH, Bayer Technology Services GmbH and the state of North- Rhine Westphalia as well 
as the European Commission in the framework of the European Regional Development 
Fund (ERDF). We also acknowledge useful discussions with Thomas Hammerschmidt, Mike 



13 



Finnis, David Pettifor and Bernd Meyer. 



* Electronic address: geor g. madsen@r ub .del 
1 W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

^ M. Mrovec, R. Groger, A. G. Bailey, D. Nguyen-Manh, C. Elsasser, and V. Vitek, Phys. Rev. 
B 75, 104119 (2007). 

^ A. Udyansky, J. von Pezold, V. N. Bugaev, M. Friak, and J. Neugebauer, Phys. Rev. B 79, 
224112 (2009). 

^ J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 

5 D. G. Pettifor and R. Podloucky, J. Phys. C.: Solid State Phys. 19, 315 (1986). 

6 A. P. Sutton, M. W. Finnis, D. G. Pettifor, and Y. Ohta, J. Phys. C.: Solid State Phys. 21, 35 
(1988). 

^ J. Harris, Phys. Rev. B 31, 1770 (1985). 

^ W. M. C. Foulkes and R. Haydock, Phys. Rev. B 39, 12520 (1989). 
9 M. W. Finnis, J. Phys.-Condes. Matter 2, 331 (1990). 

M. Elstner, D. Porezag, G. Jungnickel, J. Eisner, M. Haugk, T. Frauenheim, S. Suhai, and G. 

Seifert, Phys. Rev. B 58, 7260 (1998). 
1^ M. W. Finnis, A. T. Paxton, M. Methfessel, and M. van Schilfgaarde, Phys. Rev. Lett. 81, 5149 

(1998). 

12 D. Spanjaard and M. C. Desjonqueres, Phys. Rev. B 30, 4822 (1984). 

1^ D. Porezag, T. Frauenheim, T. Kohler, G. Seifert, and R. Kaschner, Phys. Rev. B 51, 12947 
(1995). 

1^ M. J. Mehl and D. A. Papaconstantopoulos, Phys. Rev. B 54, 4519 (1996). 

15 D. Nguyen-Manh, D. G. Pettifor, and V. Vitek, Phys. Rev. Lett. 85, 4136 (2000). 

J. Junquera, O. Paz, D. Sanchez-Portal, and E. Artacho, Phys. Rev. B 64, 235111 (2001). 
1^ A. H. Larsen, M. Vanin, J. J. Mortensen, K. S. Thygesen, and K. W. Jacobsen, Phys. Rev. B 

80, 195112 (2009). 

J. Enkovaara, C. Rostgaard, J. J. Mortensen, J. Chen, M. Dulak, L. Ferrighi, J. Gavnholt, 
C. Glinsvad, V. Haikola, H. A. Hansen, H. H. Kristoffersen, M. Kuisma, A. H. Larsen, L. 
Lehtovaara, M. Ljungberg, O. Lopez-Acevedo, P. G. Moses, J. Ojanen, T. Olsen, V. Petzold, 



14 



N. A. Romero, J. Stausholm-M0ller, M. Strange, G. A. Tritsaris, M. Vanin, M. Walter, B. 
Hammer, H. Hakkinen, G. K. H. Madsen, R. M. Nieminen, J. K. N0rskov, M. Puska, T. T. 
Rantala, J. Schi0tz, K. S. Thygesen, and K. W. Jacobsen, J. Phys.: Condens. Matter 22, 253202 
(2010). 

E. Artacho, D. Sanchez-Portal, P. Ordejon, A. Garcia, and J. Soler, Phys. Status Solidi B-Basic 
Res. 215, 809 (1999). 

20 J. J. Mortensen, L. B. Hansen, and K. W. Jacobsen, Phys. Rev. B 71, 035109 (2005). 

21 E. Francisco, L. Seijo, and L. Pueyo, J. Solid State Chem. 63, 391 (1986). 

22 D. Sanchez-Portal, E. Artacho, and J. M. Soler, J. Phys.: Condens. Matter 8, 3859 (1996). 

23 S. Kostlmeier, C. Elsasser, and B. Meyer, Ultramicroscopy 80, 145 (1999). 

24 N. Bornsen, B. Meyer, O. Grother, and M. Fahnle, J. Phys.: Condens. Matter 11, L287 (1999). 
2^ A. Urban, M. Reese, M. Mrovec, C. Elsasser, and B. Meyer, In preparation (2011). 

2^ D. G. Pettifor, Bonding and Structure of Molecules and Solids (Clarendon Press, Oxford, UK, 
1995). 

27 G. Liu, D. Nguyen-Manh, B.-G. Liu, and D. G. Pettifor, Phys. Rev. B 71, 174115 (2005). 

28 A. T. Paxton and M. W. Finnis, Phys. Rev. B 77, 024428 (2008). 

29 W. Zhong, G. Overney, and D. Tomanek, Phys. Rev. B 47, 95 (1993). 

30 D. Ycsillctcn, M. Nastar, T. A. Arias, A. T. Paxton, and S. Yip, Phys. Rev. Lett. 81, 2998 
(1998). 

31 F. Ducastelle and F. Cyrot-Lackmann, J. Phys. Chem. Solids 31, 1295 (1970). 

32 G. Allen and M. Lannoo, J. Phys. Chem. Solids 37, 699 (1976). 

33 M. W. Finnis and J. E. Sinclair, Philos. Mag. A 50, 45 (1984). 

34 M. W. Finnis, A. B. Walker, and P. Gumbsch, J. Phys.-Condcs. Matter 10, 7983 (1998). 

35 M. Mrovec, D. Nguyen-Manh, D. G. Pettifor, and V. Vitek, Phys. Rev. B 69, 094115 (2004). 

36 D. Nguyen-Manh, V. Vitek, and A. P. Horsfield, Prog. Mat. Sci. 52, 255 (2007). 

37 P.-O. Lowdin, Adv. Phys. 5, 1 (1956). 

38 B. C. Carlson and J. M. Keller, Phys. Rev. 105, 102 (1957). 

39 W. M. C. Foulkes, Phys. Rev. B 48, 14216 (1993). 

40 O. K. Andersen, Solid State Communications 13, 133 (1973). 

41 C. Kohler, G. Seifert, and T. Frauenheim, Chem. Phys. 309, 23 (2005). 

42 D. Nguyen-Manh and S. L. Dudarev, Phys. Rev. B 80, 104440 (2009). 



15 

A. T. Paxton and C. Elsasser, Phys. Rev. B 82, 235125 (2010). 

A. K. Sinha, Prog, in Mat. Sci. 15(2), 79 (1972). 

R. Drautz and D. G. Pettifor, Phys. Rev. B 74, 174117 (2006). 

D. G. Pettifor, Pliys. Rev. Lett. 63, 2480 (1989). 

M. W. Finnis, Prog. Mat. Sci. 52, 133 (2007). 



16 



Vq (A^/atom) Eq (eV/atom) Bq (GPa) c/a 



NM-FCC 

DFT 10.38 -7.890 275.59 

TB 10.38 -7.926 295.42 
NM-A15 

DFT 10.59 -7.729 271.23 

TB 10.52 -7.767 287.39 
FM-A15 

DFT 11.72 -7.978 155.05 

TB 11.90 -7.981 141.92 
NM-x 

DFT 10.55 -7.840 273.20 

TB 10.53 -7.790 271.24 
FM-BCC 

DFT 11.51 -8.064 174.38 

TB 11.58 -8.067 138.29 
AFM-FCC 

DFT 10.79 -7.946 186.42 

TB 10.74 -7.942 177.01 
NM-HCP 

DFT 10.31 -7.968 282.44 1.579 

TB 10.35 -7.966 294.54 1.570 
NM-cr 

DFT 10.55 -7.786 275.60 0.522 

TB 10.51 -7.796 267.23 0.532 



TABLE II: Equilibrium lattice constants, phase stabilities with respect to the non-magnetic free 
atom, bulk moduli and optimal c/a ratios for the studied iron compounds. 
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FE (eV) FM-BCC NM-FCC UML 

DFT 2.08 2.01 1.93 

TB (n = 0.50) 1.91 1.70 1.58 
TB (n = 0.55) 2.05 1.92 1.77 

TABLE III: Formation energies (FE) of vacancies in the FM-BCC and NM-FCC structures and of 
an unsupported monolayer of FCC-(lll) iron. The n = 0.55 model the prefactors are reoptimized 
compared to Table U giving drep — 1088 eV and <iemp — 

3.18 eV. 



